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ABSTRACT 

We report the detection of four images in the recently discovered lensed QSO 
RX J0911. 4+0551. With a maximum angular separation of 3.1", it is the quadruply 
imaged QSO with the widest known angular separation. Raw and deconvolved data 
reveal an elongated lens galaxy. The observed reddening in at least two of the four 
QSO images suggests differential extinction by this lensing galaxy. We show that both 
an ellipticity of the galaxy (e m i n = 0.075) and an external shear (7 m i n = 0.15) from 
a nearby mass has to be included in the lensing potential in order to reproduce the 
complex geometry observed in RX J0911. 4+0551. A possible galaxy cluster is detected 
about 38" from RX J091 1.4+0551 and could contribute to the X-ray emission observed 
by ROSAT in this field. The color of these galaxies indicates a plausible redshift in the 
range of 0.6-0.8. 

Subject headings: galaxies: clusters: general — gravitational lensing — quasars: 
individual (RX J0911.4+0551) 

1. Introduction 

RX J0911.4+0551, an AGN candidate selected from the ROSAT All-Sky Survey (RASS) 
(Bade et al. 1995, Hagen et al. 1995), has recently been classified by Bade et al. (1997; hereafter 
B97) as a new multiply imaged QSO. B97 show that it consists of at least three objects: two 
barely resolved components and a third fainter one located 3.1" away from the other two. They 
also show that the spectrum of this third fainter component is similar to the combined spectrum 
of the two bright components. The lensed source is a radio quiet QSO at z = 2.8. Since RASS 
detections of distant radio quiet QSOs are rare, B97 pointed out that the observed X-ray flux 
might originate from a galaxy cluster at z > 0.5 within the ROSAT error box. We present here 
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new optical and near-IR high-resolution images of RX J0911. 4+0551 obtained with the 2.56m 
Nordic Optical Telescope (NOT) and the ESO 3.5m New Technology Telescope (NTT). Careful 
deconvolution of the data allows us to clearly resolve the object into four QSO components and a 
lensing galaxy. In addition, a candidate galaxy cluster is detected in the vicinity of the four QSO 
images. We estimate its redshift from the photometric analysis of its member galaxies. 

2. Observations and reductions 

We first observed RX J0911.4+0551 in the i-C-band with IRAC 2b on the ESO/MPI 2.2m 
telescope on November 12, 1997. In spite of the poor seeing conditions (<~ 1.3"), preliminary 
deconvolution of the data made it possible to suspect the quadruple nature of this object. Much 
better optical observations were obtained at the NOT (La Palma, Canary Islands, Spain). Three 
300s exposures through the / filter, with a seeing of ~ Of! 8 were obtained with ALFOSC under 
photometric conditions on November 16, 1997. Under non-photometric, but excellent seeing 
conditions (~ Of! 5 — Of! 8), three 300s /-band exposures, three 300s U-band and five 600s [/-band 
exposures were taken with HIRAC on the night of December 3, 1997. The pixel scales for HIRAC 
and ALFOSC are 0'.'1071 and / . / 186, respectively. RX J0911.4+0551 was also the first gravitational 
lens to be observed with the new wide field near-IR instrument SOFI, mounted on the ESO 
3.5m NTT. Excellent K and J images were taken on December 15, 1997, and January 19, 1998 
respectively. The 1024 x 1024 Rockwell detector was used with a pixel scale of 0'.'144. 

The optical data were bias subtracted and flat-field corrected using sky- flats. Fringe-correction 
was also applied to the I-band data from ALFOSC. Sky subtraction was carried out by fitting 
low-order polynomial surfaces to selected areas of the frames. Cosmic ray removal was finally 
performed on the data. The infrared data were processed as explained in Courbin, Lidman & 
Magain (1998), but in a much more efficient way for SOFI than for IRAC 2b data, since the array 
used with the former instrument is cosmetically superior to the array used with the latter. 

3. Deconvolution of RX J0911.4+0551 

All images were deconvolved using the new deconvolution algorithm by Magain, Courbin 
& Sohy (1998; hereafter MCS). The sampling of the images was improved in the deconvolution 
process, i.e. the adopted pixel size in the deconvolved image is half the pixel size of the original 
frames. The final resolution adopted in each band was chosen according to the signal-to-noise 
(S/N) ratio of the data, the final resolution improving with the S/N. 

Our NOT HIRAC data were deconvolved down to the best resolution achievable with the 
adopted sampling (2 pixels FWHM, i.e. 0f!l), whereas the NOT/ALFOSC frames were deconvolved 
to a resolution of 3 pixels FWHM, i.e. Of! 29. Although of very good quality, the near IR-data have 
a lower S/N than the optical data. The resolution was therefore limited to 5 pixels FWHM in 
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both J and K ((X'36) in order to avoid noise enhancement. 

The MCS algorithm can be used to deconvolve simultaneously a stack of individual images 
or to deconvolve a single stacked image. The results from the simultaneous deconvolutions are 
displayed in Fig. |l[ The four QSO images are labeled Al, A2, A3, and B. The quality of the results 
was checked from the residual maps, as explained in Courbin et al. (1998). The deconvolution 
procedure decomposes the images into a number of Gaussian point sources, for which the program 
returns the positions and intensities, plus a deconvolved numerical background. In the present 
case, the deconvolution was also performed with an analytical De Vaucouleurs, and exponential 
disk galaxy profile at the position of the lensing galaxy, in order to better describe its morphology. 

3.1. Results 

Table 1 lists the flux of each QSO component, relative to Al, as derived from the simultaneous 
deconvolutions. Although the HIRAC U, V and / band data were taken during non photometric 
conditions, they can still be used to determine the relative fluxes between the four images of the 
QSO. The errors in the relative fluxes are determined from the simultaneous deconvolutions and 
represent the 1-a standard deviation in the peak intensities. 

When flux calibration was possible, magnitudes were calculated and the results are displayed 
in Table 2, which also contains the positions of the four QSO components relative to Al. The 
astrometric errors are derived by comparing the positions of the components in the different bands. 

The deconvolved numerical background is used to determine the galaxy position from the 
first order moment of the light distribution. A reasonable estimate of the error on the lens 
position was derived from moment measurements through apertures of varying size and on 
several images obtained by running deconvolutions with different initial conditions. /, J, and 
K magnitudes for the galaxy were also estimated from the deconvolved background image by 
aperture photometry (~ 1.2" aperture in diameter). The numerical galaxy is elongated in all 
three bands and the position angle of its major axis is Qq ~ 140 ±5°. In the near-IR, the galaxy 
looks like it is composed of a bright sharp nucleus plus a diffuse elongated disk. However we can 
not exclude that the observed elongation is due to an unresolved blend of two or more intervening 
objects. None of our two analytical profiles fit perfectly the galaxy. The De Vaucouleurs profile 
(e = 1 - b/a = 0.31 ± 0.07, 9 G = 130 ± 10°) fits slightly better than the exponential disk light 
distribution, but still produces residual maps with values as large as 1.5-2 per pixel, compared 
with a x 2 °f 2.5-5 per pixel for the exponential disk profile. The ellipticity and the position angle 
derived this way are very uncertain due to the low S/N of the data. Much deeper observations will 
be required to perform precise surface photometry of the lens(es) and to draw a definite conclusion 
about its (their) morphology. 
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3.2. Field photometry 

In order to detect any intervening galaxy cluster which might be involved in the overall 
lensing potential and contributing to the X-ray emission observed by ROSAT, we performed J, J, 
and K band photometry on all the galaxies in a 2.5' field around the lensed QSO. A composite 
color image was also constructed from the frames taken through these 3 filters, in order to directly 
visualize any group of galaxies with similar colors, and therefore likely to be at the same redshift. 
The color composite is presented in Fig. |2[ 

Aperture photometry was carried out using the SExtractor package (LINUX version 1.2b5, 
Bertin & Arnouts, 1996). The faintest objects were selected to have at least 5 adjacent pixels 
above 1.2a s k y , leading to the limiting magnitudes 23.8, 21.6, and 20.0 mag/arcsec 2 in the /, J, 
and K bands, respectively. The faintest extended object measured in the different bands had 
magnitudes 23.0, 22.0, and 20.3 in /, J, and K, respectively. The color-magnitude diagram of the 
field was constructed from the I and K band data which give the widest wavelength range possible 
with our photometric data. Since the seeing was different in the two bands, particular attention 
was paid to the choice of the isophotal apertures fitted to the galaxies. They were chosen to be 
as large as possible, still avoiding too much contamination from the sky noise, as an oversized 
isophotal aperture would introduce. The color-magnitude diagram of the galaxies in the 2.5' field 
is displayed in Fig. |3] . Stars are not included in this plot. Note that, because of their proximity, 
the magnitudes of the two blended members of the cluster candidate (center of the circle in Fig. |2|) 
might be underestimated by as much as 0.3-0.4 magnitudes in both I and K. 



4. Models 

In a first attempt to model the system we chose an elliptical potential of the form 

V> = V^ e ), (1) 
2 x' 2 y' 2 

Te ~ (1 + e) 2 + (1-e) 2 ' () 

where the coordinates x' and y' are measured along the principal axes of the galaxy, whose position 
angle 9q is a free parameter. For small ellipticities e, this potential is a good approximation for 
elliptical mass distributions (Kassiola & Kovner 1993). Additionally, an external shear 7 with 
direction <f> is included. Even without assuming an explicit potential if), we can determine the 
minimal 7 and e needed to reproduce the observed image positions, by applying the methods 
given by Witt & Mao (1997). Their methods eliminate the unknown parameters of the model to 
find constraints for the ellipticity of the galaxy and the external shear. For the shear we predict 
a minimum of 7min = 

0.15 ± 0.07, while 

e min — 0.075 i 0.034 (both with lcr errors). Therefore, 
neither ellipticity nor shear should be omitted from the modeling. To keep models simple, we used 
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a pseudo-isothermal potential corresponding to an apparent deflection angle of 

_ QQ f x y \ 

Q r c + ^f+^ V(l + e) 2 ' (l-e)V' U 

which degenerates to a singular isothermal sphere for vanishing core-radius r c and ellipticity e. 

The model's parameters are determined by minimizing the x 2 given by the observed and 
predicted positions of the images and the galaxy. The best model leads to x 2 = 0-65 and has 
parameters as given in Table ||. This value is reasonable for a model with one degree of freedom 
(r c can not be counted as a free parameter here, because without the restriction of r c > it would 
become negative in the fit). The positions of the images and the galaxy in Table || are parameters 
of the model and have to be compared with the observed positions in Table ^. 

As shown above, no elliptical potential without a shear, and no spherical potential with shear 
can reproduce the observational data. For our pseudo-isothermal model, this results in large 
X 2 values of 50.5 and 63.7, respectively. These values are much too high for the two degrees of 
freedom. 



5. Discussion 

Thanks to our new high-resolution imaging data, the QSO RX J0911. 4+0551 is resolved into 
four images. In addition, deconvolution with the new MCS algorithm reveals the lensing galaxy, 
clearly confirming the lensed nature of this system. The image deconvolution provides precise 
photometry and astrometry for all the components of the system. 

Reddening in components A2 and A3 relative to Al is observed from our U, V, and / frames 
that were taken within three hours on the same night. The absence of reddening in component 
B and the difference in reddening between components A2 and A3 suggest extinction by the 
deflecting galaxy. Note that although our near-IR data were obtained from 15 days to 6 weeks 
after the optical images, they appear to be consistent with the optical fluxes measured for the 
QSO images, i.e. flux ratios increase continuously with wavelength, from U to K, indicating 
extinction by the lensing galaxy. 

We have discovered a good galaxy cluster candidate in the SW vicinity of RX J0911. 4+0551 
from our field photometry in the /, J, and K bands. Comparison of our color-magnitude 
diagram with that of a blank field (e.g., Moustakas et al. 1997) shows that the galaxies around 
RX J0911. 4+0551 are redder than field-galaxies at an equivalent apparent magnitude. In addition, 
the brightest galaxies in Fig. || lie on a red sequence at I — K ~ 3.3, typical for the early type 
members of a distant galaxy cluster. The two dashed lines indicate our ±0.4 color error bars at 
K ~ 19 around I — K ~ 3.3. Most of these galaxies are grouped in the region around a double 
elliptical at a distance of ~ 38" and a position angle of ~ 204° relative to Al. This can also 
be seen in Fig. || which shows a group of red galaxies with similar colors centered on the double 
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elliptical (in the center of the circle). 

Consequently, there is considerable evidence for at least one galaxy cluster in the field. The 
redshift of our best candidate cluster (the one circled in Fig. |2|) can be estimated from the I and 
K band photometry. We have compared the -fT-band magnitudes of the brightest cluster galaxies 
with the empirical K magnitude vs. redshift relation found by Aragon-Salamanca et al. (1998). 
We find that our cluster candidate, with its brightest K magnitude of about ~ 17.0, should have 
a redshift of z ~ 0.7. A similar comparison has been done in the /-band without taking into 
account galaxy morphology. We compare the mean I magnitude of the cluster members with the 
ones found by Koo et al. (1996) for galaxies with known redshifts in the Hubble Deep Field and 
obtain a cluster redshift between 0.6 and 0.9. Finally, comparison of the I — K color of the galaxy 
sequence with data and models from Kodama et al. (1998) confirm the redshift estimate of 0.6-0.8. 

In order to calculate physical quantities from the model parameters found in section 4, we 
assume a simple model for the cluster which may be responsible for the external shear. For an 
isothermal potential, the true shear and convergence are of the same order of magnitude. As 
the convergence is not explicitly included in the model, the deduced shear is a reduced shear 
leading to an absolute convergence of k = 7/(1 + 7) = 0.241. For a cluster redshift of = 0.7 
and with cosmological parameters $1 = 1, A = this corresponds to a velocity dispersion of about 
1100 km s -1 if the cluster is positioned at an angular distance of 40" . See Gorenstein, Falco Sz 
Shapiro (1988) for a discussion of the degeneracy preventing a direct determination of k. From the 
direction of the shear <f>, (see Table |||) we can predict the position angle of the cluster as seen from 
the QSO to be 12° or 192°. The latter value agrees well with the position of our cluster candidate 
SW of the QSO images. Note also the good agreement between the position angle 0q derived from 
the observed light distribution, and the predicted position angle corresponding to our best fitting 
model of the lensing potential. Interestingly, this is in good agreement with Keeton, Kochanek & 
Falco (1998) who find that projected mass distributions are generally aligned with the projected 
light distributions to less than 10°. 

The color of the main lensing galaxy is very similar to that of the cluster members, suggesting 
that it might be a member of the cluster. Using the same model for the cluster as above, assuming 
the galaxy at the same redshift as the cluster, and neglecting the small ellipticity of e < 0.05, the 
velocity dispersion of the lensing galaxy can be predicted from the calculated deflection angle ao 
to be of the order of 240 km s . Since the galaxy profile is sharp towards the nucleus in K, we 
cannot rule out the possibility of a fifth central image of the source, as predicted for non-singular 
lens models. Near-IR spectroscopy is needed to get a redshift determination of the lens and to 
show whether it is blended or not with a fifth image of the (QSO) source. 

Some 10" SW from the lens, we detect a small group of even redder objects. These red 
galaxies can be seen in Fig. [2] a few arcseconds to the left and to the right of the cross. They might 
be part of a second galaxy-group at a higher redshift, and with a position in better agreement 
with the X-ray position mentioned by B97. However, since the measured X-ray signal is near the 
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detection limit, and the 1-a positional uncertainty is at least 20" , the X-ray emission is compatible 
with both the QSO and these galaxy groups in the field. Furthermore, this second group, at 
z > 0.7, would most likely be too faint in the X-ray domain to be detected in the RASS. In fact, 
even our lower redshift cluster candidate would need to have an X-ray luminosity of the order of 
Lo.i-2.4keV ~ 7.10 44 ergs _1 (assuming a 6 keV thermal spectrum, Ho = SOkms" 1 Mpc -1 , qo = 0.5), 
in order to be detected with 0.02 ctss -1 by ROSAT. This is very bright but not unrealistic for 
high redshift galaxy clusters (e.g., MS 1054-03, Donahue, Gioia, Luppino et al. 1997). 

RX J0911. 4+0551 is a new quadruply imaged QSO with an unusual image configuration. The 
lens configuration is complex, composed of one main lensing galaxy plus external shear possibly 
caused by a galaxy cluster at redshift between 0.6 and 0.8 and another possible group at z > 0.7. 
Multi-object spectroscopy is needed in order to confirm our cluster candidate/s and derive its/their 
redshift and velocity dispersion. In addition, weak lensing analysis of background galaxies might 
prove useful to map the overall lensing potential involved in this complex system. 
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photometric system, T.J. Broadhurst, C.R. Keeton, G. Rauw for helpful discussions and R. Sip for 
help in improving the figures. IB, FC and R0 are supported in part by contract ARC94/99-178 
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Table 1: Flux ratios for all four QSO components, relative to Al 



A2 A3 B 

~F^ 0.965 ±0.013 0.544 ± 0.025 0.458 ± 0.004 

Fj 0.885 ±0.003 0.496 ± 0.005 0.412 ± 0.005 

Fj 0.680 ±0.013 0.398 ±0.002 0.420 ± 0.003 

F v 0.587 ±0.009 0.334 ± 0.004 0.413 ± 0.006 

F v 0.590 ±0.013 0.285 ±0.007 0.393 ± 0.004 



Note. — Results obtained from our non-photometric data. All measurements are given along with their 1-a errors. 
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Fig. 1. — High resolution images of RX J0911.4±0551 are displayed in the left panels with, in the 
right panels, their corresponding deconvolved versions with improved resolution and sampling. 
In each band, the individual data frames are deconvolved simultaneously. [Top:] Stack of 3 
NOT±HIRAC I-band image. The total exposure time is 900s. [Middle:] Stack of 9 NTT±SOFI 
J-band frames. The total exposure time is 1080s. [Bottom:] Stack of 19 NTT±SOFI K-band with 
a total exposure time of 2400s. For all the three bands the object is clearly resolved into four QSO 
images, labeled Al, A2, A3, and B, plus the elongated lensing galaxy. The field of the optical and 
near IR data are respectively 7" and 9" on a side. North is to the top and East to the left in all 
frames. 
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Fig. 2. — Composite image of a 2' field around RX J0911. 4+0551. The frame has been obtained 
by combining our /, J and if -band data. North is up and East to the left. Note the group of red 
galaxies with similar colors, about 38" SW of the quadruple lens (circle) and the group of even 
redder galaxies 10" SW of the lens (cross). 




Fig. 3. — Color magnitude diagram of the 2.5' field around RX J0911. 4+0551. The lensing galaxy 
and many of the cluster members form the overdensity around / — K ~ 3.3. The lensing galaxy 
in RX J0911. 4+0551 is marked by a star and the two blended galaxies in the center of the cluster 
candidate are plotted as triangles. Stars are not plotted in this diagram. 
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T able 2: Photometric and astrometric properties of RX J0911. 4+0551 and the lensing galaxy. 





Al 


A2 


A3 


B 


G 


K 


17.137 ±0.020 


17.308 ± 0.035 


17.933 ± 0.074 


18.119 ±0.010 


17.88 ±0.12 


J 


17.794 ± 0.030 


17.927 ± 0.029 


18.556 ±0.038 


18.757 ±0.020 


19.85 ±0.12 


I 


18.310 ±0.084 


18.728 ± 0.096 


19.305 ±0.133 


19.197 ±0.083 


21.36 ±0.13 


x(") 


0.000 ± 0.004 


-0.259 ± 0.007 


±0.013 ±0.008 


±2.935 ± 0.002 


±0.709 ± 0.026 


y(") 


0.000 ± 0.008 


±0.402 ± 0.006 


±0.946 ± 0.008 


±0.785 ± 0.003 


±0.507 ±0.046 



Note. — The astrometry is given relative to component Al with x and y coordinates defined positive to the North 
and West of Al. All measurements are given along with their 1-cr errors. The 1-a errors on the photometric zero 
points are 0.02 in all bands. 



Table 3: Parameters of the best-fit model. 



«o r c e 9 G 
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0s 


l'/1061 0.0476 140?5 


0.3176 


102?3 ±0"7078 


±1'/1459 






±0'/5440 


±0'/5452 



image 


X 


y 


amplification 


Ai 


±0'/0000 


-0'/0004 


-4.45 


A 2 


-0'.'2590 


±0'/4019 


±8.59 


A 3 


±0'/0130 


±0'/9456 


-3.70 


B 


±2'/9350 


±0'/7850 


±1.79 



Note. — 0o and Gs are the position of the galaxy and the source, respectively. The core radius r c was a free 
parameter constrained to r c > in the modeling. The best fit lead to r c — 0. Position angles 9q and <f> are measured 
from N to E. Negative amplification is due to odd parity of the image. 
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